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We have searched for the lepton-number violating processes B"*" — ^ h~ t'^t'^ with h~ = K~ /tt~ 
and = e^ I jjL^ , using a sample of 471 ± 3 million BB events collected with the BABAR detector 
at the PEP-II e^e~ collider at the SLAC National Accelerator Laboratory. We find no evidence 
for these decays and place 90% confidence level upper limits on their branching fractions B [B'^ — ^ 
7r-e+e+) < 2.3 x 10-^ B (B+ -> A'-e+e+) < 3.0 x 10"*, B {B+ Tr-fi+fi+) < 10.7 x 10"*, and 
B {B+ K-fi+fi+) < 6.7 X 10"*. 

PACS numbers: 13.20.He,13.15.+g,14.60.St 



In the Standard Model (SM) , Icpton number L is con- 
served in low-energy collisions and decays [ij and the lep- 
ton flavor numbers for the three lepton families are con- 
served if neutrinos are massless. The observation of neu- 
trino oscillations 0] indicates that neutrinos have mass. 
If the neutrinos are of the Major ana type Q, the neutrino 
and antineutrino are the same particle and processes that 
involve lepton-number violation become possible. The 
lepton number must change by two units [AL = 2) in 
this case and the most sensitive searches have so far in- 
volved neutrinoless nuclear double beta decays Ovfip [1]. 
The nuclear environment complicates the extraction of 
the neutrino mass scale. Processes involving meson de- 
cays have been proposed as an alternative that can also 
look for lepton-number violation with muons or r Icp- 
tons. 

An example of a decay involving mesons is -B^ — > 
where = e+ or yu+ and /i" is a meson with a 
mass smaller than the B meson. A possible mechanism 
for this process involving the production and subsequent 
decay of a Majorana neutrino is illustrated in Fig. [TJ 
which is topologically similar to the t-channel Feynman 
diagram in Oi^/3/3 decays. If the Majorana neutrino mass 
lies between the h meson and the B meson masses, res- 
onance production could result in an enhanced peak in 



the invariant mass spectrum of the hadron and one of the 
leptons 
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FIG. 1: An example diagram of the AL = 2 process B'^ — > 
h~ t'^ via s-channel Majorana neutrino Vm production and 
decay. 



The experimental approach in searches for these 
lepton-number violating processes is very similar to the 
approach for B — >■ K'^*^t^(.~ and we use many of the tech- 
niques reported in Refs. [1,[3| to search for the four modes 
£?+ ji^ i^L^ and h^e^e~^, where h~ ~ or 

TT" (charge-conjugate modes are implied throughout this 
paper) . Previous searches for these decays have produced 
90% confidence level (C.L.) upper limits on the branch- 
ing fractions in the range (1.0 - 8.3) x lO"*^ The 
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LHCb collaboration recently reported 95% C.L. upper 
limits on the branching fractions B{B^ — K~ fi^ ^~^) < 
5.4 X 10-« and B{B+ Tr-fi+fi+) < 1.3 x IQ-^ 
The Belle collaboration place 90% C.L. upper limits on 
the branching fractions B{B^ — > D^i^i^) in the range 
(1.1 - 2.6) X 10-6 [13. 

Our search uses a data sample of 471 ± 3 million BB 
pairs collected at the T{AS) resonance with the BaBAR 
detector at the PEP-II asymmetric-energy e+e" collider 
at the SLAG National Accelerator Laboratory. The 
e~^e~ center-of-mass (CM) energy is y/s = 10.58 GeV, 
corresponding to the mass of the T{AS) resonance (on- 
resonance data). In addition, 43.9 fb~^ of data collected 
40 MeV below the T(4S') resonance (off-resonance data) 
are used for background studies. We assume equal pro- 
duction rates of B^B~ and i?°B° mesons lla . The 
BaBAR detector is described in detail in Ref. [ll| . 

Monte Carlo (MC) simulation is used to evaluate the 
background contamination and selection efficiencies. The 
simulated backgrounds are also used to cross-check the 
selection procedure and for studies of systematic effects. 
The signal channels are simulated by the EvtGen [l3| 
package using a three-body phase space model. We also 
generate light quark qq continuum events (e'^e^ qq, 
q = u,d,s,c), di-muon, Bhabha elastic e"'"e~ scattering, 
BB background and two-photon events [31 . Final-state 
radiation is provided by Photos 15 1. The detector re- 
sponse is simulated with GEANT4 [l6|, and all simulated 
events are reconstructed in the same manner as data. 

We select events that have at least four charged tracks, 
the ratio of the second to zeroth Fox- Wolfram mo- 
ments 17[ less than 0.5, and two same-sign charged lep- 



tons each with momentum greater than 0.3 GeV/c in the 
laboratory frame. The total transverse vector momen- 
tum of an event calculated in the laboratory frame must 
be less than 4GeV/c; the distribution of this quantity 
peaks at 0.2GeV/c for signal events. The two leptons 
are constrained to come from a single vertex and an in- 
variant mass m£+£+ < 5.0GcV/c^ is required, to main- 
tain compatibility with Ref. Q. Electrons and positrons 
from photon conversions are removed, where photon con- 
version is indicated by electron-positron pairs with an 
invariant mass less than 0.03 GeV/ and a production 
vertex more than 2 cm from the beam axis. 

The charged pions and kaons are identified by mea- 
surements of their energy loss in the tracking detectors, 
the number of photons recorded by the ring-imaging 
Cherenkov detector and the corresponding Chercnkov 
angle. These measurements are combined with informa- 
tion from the electromagnetic calorimeter and the instru- 
mented magnetic-flux return detector to identify elec- 
trons and muons 11|. 

The four-momenta of the electrons and positrons are 
corrected for Bremsstrahlung emission by searching for 
compatible photons. Using measurements made in the 
laboratory frame, the photon and electron four-momenta 



arc combined if the photon energy Ej is greater than 
0.05 GeV, the shape of the energy deposit in the electro- 
magnetic calorimeter is compatible with a photon shower, 
and the difference in polar angle between the photon and 
electron, measured at the point of closest approach to 
the beam spot, is less than 0.035 rad. In addition, the 
azimuthal angles (j) of the photon the lepton (/)£, and 
the calorimeter deposit associated with the lepton (pc, all 
measured at the primary vertex, must be compatible with 
0^ — 0.05 < < (j)c for electrons and (j)c < <j)j < 4>e + 0-05 
for positrons. 

The two leptons and the hadron track arc combined 
to form a signal B candidate. The B candidate is re- 
jected if the invariant mass of the two leptons is in the 
range 2.85 < mi+i+ < 3.15 GeV/ or 3.59 < m(+i+ < 
3.77 GeV/ c^. Although a peaking background in the 
J/ip or il){2S) mass regions is not expected, these cri- 
teria maintain consistency with Ref. [6|. For the mode 
B^ — > TT~ n'^ fx'^ , the invariant mass of each muon and 
the hadron must be outside the region 3.05 < m^+j^- < 
3.13 GeV/ c^. This rejects events where a muon from a 
J/ip decay is misidcntified as a pion. The probability to 
misidentify a pion as a muon is of the order 2% and to 
misidentify as an electron less than 0.1%. 

We measure the kinematic variables rngs = 
\/s/A — p*g and AE ~ Eg — -^5/2, where p% and Eg are 
the B momentum and energy in the T(45') CM frame, 
and -y/s is the total CM energy. For signal events, the 
mEs distribution peaks at the B meson mass with a 
resolution of about 2.5MeV/c^ and the AE distribution 
peaks near zero with a resolution of about 20 MeV, in- 
dicating that the candidate system of particles has to- 
tal energy consistent with the beam energy in the CM 
frame. The signal B candidate is required to be in 
the kinematic region 5.200 < tties < 5.289 GeV/c^ and 
-0.10 <AE< 0.05 GeV. 

The main backgrounds arise from light quark qq con- 
tinuum events and BB backgrounds formed from ran- 
dom combinations of leptons from semileptonic B and 
D decays. These are suppressed through the use of 
boosted decision tree discriminants (BDTs) [3] ■ As the 
input variable distributions for the qq continuum and the 
BB backgrounds are sufficiently different, two BDTs are 
trained, one to distinguish between signal and qq contin- 
uum, and the other between signal and BB backgrounds. 
Each BDT is trained in four regions according to lep- 
ton type (muon versus electron) and mass range {mg+(+ 
above or below the J/ip mass). The input variables con- 
sist of AE and seventeen parameters that represent the 
event shape of the decay, the distance of closest approach 
of the di-lepton system to the beam axis, the vertex prob- 
abilities of the di-lepton and B candidates, the magni- 
tudes of the thrusts of both the decay particles and the 
rest of the event, and the thrust directions with respect 
to the beam axis of the experiment. 

To construct the BDTs, we use simulated samples of 
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events for the signal and background, and we assume 
background decay rates consistent with measured val- 
ues [l9|. We compare the distributions of the data and 
the simulated background variables used as input to the 
BDTs and confirm that they are consistent. 

The output distributions of the qq and BB BDTs are 
each used to define probability distribution functions Vsig 
and Pbkg for signal and background, respectively. The 
probabilities are used to define a likelihood ratio H. as 
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We veto candidates if either V^J^ or V!!Z is less than 0.5 
or the ratio TZ less than 0.2. This retains 85% of the 
simulated signal events while rejecting more than 95% of 
the background. 

After the application of all selection criteria, some 
events will contain more than one reconstructed B can- 
didate. Less than 1% of accepted events have more than 
one signal candidate. We select the most probable signal 
candidate from among all the candidates in the event us- 
ing the likelihood ratio TZ. Averaged over all events, the 
correct candidate in simulated signal events is selected 
with greater than 98.5% accuracy. For events with more 
than one candidate, the correct candidate is selected with 
an accuracy of 67%-82% depending on the mode. The fi- 
nal event selection efficiency for simulated signal is 13%- 
48% depending on the final state. The selection efficiency 
for all modes is approximately constant to within a rel- 
ative ±10% as a function of to^+^j- between to^- and 
4.6GeV/c2. 

We extract the signal and background yields from the 
data with an unbinned maximum likelihood (ML) fit us- 
ing 
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where the likelihood C for each event candidate i is the 
sum of njT'j(xi] dj) over two categories j: the signal 
mode 5+ h~i^£^ (including the small number of 
misreconstructed signal candidates) and background as 
discussed below. For each category j, T'j{xi;dj) is the 
product of the probability density functions (PDFs) eval- 
uated for the i-th event's measured variables Xi. The 
number of events for category j is denoted by rij and TV 
is the total number of events in the sample. The quanti- 
ties dj represent the parameters describing the expected 
distributions of the measured variables for each category 
j. Each discriminating variable Xi in the likelihood func- 
tion is modeled with a PDF, where the parameters dj 
are extracted from MC simulation, off-resonance data. 



or on-resonance data with tties < 5.27GeV/c^. The two 
variables Xi used in the fit are uies and TZ. Since the lin- 
ear correlations between the two variables are found to 
be only 4%-7% for simulated signal modes and 8%-12% 
for simulated background and on-resonance data, we take 
each T'j to be the product of the PDFs for the separate 
variables. Any correlations in the variables are treated 
later as a systematic uncertainty. The three free param- 
eters in the fit are the numbers of signal and background 
events and the slope of the background tties distribution. 

MC simulations show that the qq and BB backgrounds 
have very similar distributions in mes and TZ. We there- 
fore use a single ARGUS shape [2^ to describe the toes 
combinatorial background, allowing the shape parameter 
to float in the fits. The ratio TZ for both signal and back- 
ground is fitted u sing a non-parametric kernel estimation 
KEYS algorithm [2l[. 

We parameterize the signal toes distributions using a 
Gaussian shape unique to each final state, with the mean 
and width determined from fits to the analogous final 
states in the £?+ — > J/i/'(— > £+^^)/i+ events from the 
on-resonance data. The same selection criteria as given 
above are used, with the modification that two opposite- 
sign leptons are required, the reconstructed J/ip mass 
must be in the range 2.95 to 3.15GeV/c^, toes greater 
than 5.24GeV/c2, and AE between -0.3 and 0.2 GeV. 
The signal and background toes on-resonance data dis- 
tributions are fitted with a Gaussian and an ARGUS 
function, respectively. For modes with a pion in the final 
state, we account for J/-0if+ misidentified as J/ip7T~^ by 
using the signal distribution extracted from the J/tjjK^ 
data as an additional background. For both J/ip — > e^e^ 
and J/^p — > 11'^ , the J/^p mass distribution has a 
width ~' 15McV/c^. The toes mean for all modes is 
5.2791 ± 0.0001 GcV/c2 and the width is (2.41 - 2.56) 
MeV/c2 with an error (0.02 - 0.09) MeV/c^, depending 
on the mode. The means and widths are robust against 
changes in the assumptions concerning the relative con- 
tribution of qq and BB events to the backgrounds and 
the functions used to fit the signal and background dis- 
tributions. The numbers of measured events for the four 
modes are within one standard deviation of the expected 
numbers calculated from previously measured branching 
fractions 12| . Figure [2] shows the extracted toes distri- 
butions for each mode. 

As a cross-check of the background PDFs to B^ ^■ 
h~£~^£~^ , we perform a fit to a simulated background sam- 
ple, with the same number of events as the on-resonance 
data sample, and also a fit to the off-resonance data sam- 
ple. In both cases, the number of signal events is com- 
patible with zero for all four modes. 

We test the performance of the fits to i?+ — )- h~i^£~^ 
by generating ensembles of MC datasets from both the 
PDF distributions and the fully simulated MC events. 
The mean number of signal and background events used 
in the ensembles is taken from the full default model fit to 
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FIG. 2: The tues distributions for a) B+ -s- J/ipi^ 
e+e~)K+- b) B+ J/^ iJ,+ fj.-)K+- c) B+ J/->p {^ 
e+e")7r+; and d) B+ -> J/^ (-;■ At+M")7r+. The solid (blue) 
line is the total fit, the dotted (magenta) line is the back- 
ground, the dash-dotted (green) line is the signal, the dashed 
(black) line is the misidentified J/tpK'^ events. 



FIG. 3: The rriBS distributions for a) B+ -i> K-e+e+; b) 
B+ K-jj,+fj,+ ; c) B+ n~e+e+; and d) B+ n- fi+ fi+ . 
The solid (blue) line is the total fit, the dotted (magenta) line 
is the background, the solid (green) histogram is the signal. 



the selected on-resonance data sample described above. 
We generate and fit 5000 datasets with the number of sig- 
nal and background events allowed to fluctuate according 
to a Poisson distribution. The signal yield bias in the en- 
semble of fits is between -0.30 and 0.15 events, depending 
on the mode, and this is subtracted from the yield taken 
from the data. 

The results of the ML fits to the on-rcsonance data 
arc summarized in Table ID Figure [3] shows the toes dis- 
tributions for the four modes. The signal significance is 
defined as S = \/2Aln£, where Aln£ is the change in 
log-likelihood from the maximum value to the value when 
the number of signal events is set to zero. Systematic 
errors are included in the log-likelihood distribution by 
convolving the likelihood function with a Gaussian dis- 
tribution with a variance equal to the total systematic 
error defined below. The branching fraction B is given 
by ris/ [rjN where is the signal yield, corrected for 
the fit bias, 77 is the reconstruction efficiency and N is 
the number of BB events collected. 

The systematic uncertainties on the branching frac- 
tions are summarized in Table |IT1 They arise from the 
PDF parameterization, fit biases, background yields, and 
efficiencies. The PDF uncertainties are calculated by 
varying, by their errors, the PDF parameters that are 
held fixed in the default fit, taking into account correla- 
tions. For the non-parametric kernel estimation KEYS 
algorithm, we vary the smearing parameter between 50% 
and 200% of the nominal value. The uncertainty from 
the fit bias includes the statistical uncertainty from the 
simulated experiments and half of the correction itself, 



added in quadrature. 

Two tests are used to calculate the contribution to the 
error caused by the assumption that the qq and BB back- 
grounds have similar distributions. We first vary the rel- 
ative proportions of light quark qq, cc, and BB used in 
the simulated background between 0% and 100%. The 
new simulated background TZ PDF is then used in the 
fit to the data and compared to the default fit to data. 
We also perform an ensemble of fits to MC samples con- 
sisting of one simulated signal event and the number of 
simulated background events given by the default fit to 
data. The relative proportions of light quark qq, cc, and 
BB in the simulated background are varied and a fit per- 
formed to the MC sample. The result is compared to the 
fit to the default MC sample. The error is calculated as 
half the difference between the default fit and the max- 
imum deviation seen in the ensemble of fits. All these 
errors described above are additive in nature and affect 
the significance of the branching fraction results. 

Multiplicative uncertainties include reconstruction effi- 
ciency uncertainties from tracking (0.8% per track added 
linearly for the leptons and 0.7% for the kaon or pion), 
charged Icpton particle identification (0.7% per track 
added linearly for electrons, 1.0% for muons), hadrou 
particle identification (0.2% for pions, 0.6% for kaons), 
uncertainty in the BDT response from comparison to 
charmonium control samples (2.0%), the number of BB 
pairs (0.6%), and MC signal statistics (0.2%). The total 
multiplicative branching fraction uncertainty is 3.2% or 
less for all modes. 

As shown in Table HI we observe no significant yields. 
The 90% C.L. branching fraction upper limits Bjjl are 
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TABLE I: Results for the measured B decays, showing the total events in the sample, signal yield fit bias (with error), signal 
yield (corrected for fit bias) and its statistical uncertainty, reconstruction efficiency rj, significance S (with statistical and 
systematic uncertainties included), branching fraction B, and 90% C.L. branching fraction upper limits Bul- 

Mode Events Fit Bias Yield 77 (%) S ja) B (xlO~^) Bul (x10~^) 

B+~^-K-e+e+ 123 +0.15 ± 0.09 O.Bta'r 47.8 ±0.1 0.4 0.271^ 2+0.1 2.3 

B+^K-e+e+ 42 -0.30 + 0.15 0.7+}:2 30.9 + 0.1 0.5 0.49ti:g+0.1 3.0 

B+ TT- fj,+ fi+ 228 -0.01 + 0.05 O.Otl l 13.1 +0.1 0.0 0.031^:2+0.6 10.7 

B+~^K-fi+fi+ 209 +0.02 + 0.04 O-Sttl 23.0 + 0.1 0.2 0.451^:7+0.4 6.7 



TABLE IL Summary of branching fraction B systematic un- 
certainties for the four decays. 



Systematic tt e'^e'^ K n fi'^ j-i'^ K' 



Additive uncertainties (candidates) 



PDF variation 


0.01 


0.01 


0.01 


0.09 


KEYS PDFs 


0.30 


0.05 


0.23 


0.02 


Fit bias 


0.09 


0.15 


0.05 


0.04 


Backgrounds 


0.05 


0.07 


0.25 


0.35 


Total 


0.32 


0.17 


0.34 


0.35 



Multiplicative uncertainties (%) 



Lepton tracking 


1.6 


1.6 


1.6 


1.6 


Hadron tracking 


1.4 


1.4 


1.4 


1.4 


Lepton ID 


0.7 


0.7 


1.0 


1.0 


Hadron ID 


0.2 


0.6 


0.2 


0.6 


BDT 


2.0 


2.0 


2.0 


2.0 


BB pairs 


0.6 


0.6 


0.6 


0.6 


MC statistics 


0.2 


0.2 


0.2 


0.2 


Total 


3.1 


3.1 


3.2 


3.2 



Branching 


fraction B uncertainties (xlO 




Additive 


0.14 


0.12 


0.56 


0.34 


Multiplicative 


0.01 


0.02 


0.01 


0.01 


Total 


0.14 


0.12 


0.56 


0.35 



determined by integrating the total likelihood distribu- 
tion (taking into account statistical and systematic er- 
rors) as a function of the branching fraction from to 
Bul, such that ^^""^ CdB = 0.9 x CdB. The upper 
limits are dominated by the statistical error. 

FigurelHshows the 90% C.L. upper limit on the branch- 
ing fractions as a function of the mass mg+f^- for the four 
modes. The limit is recalculated in bins of O.lGeV/c^ 
with the assumption that all the fitted signal events are 
contained in that bin. The change in shape is due to 
the variation of the reconstruction efficiency as a func- 
tion of the mass. If the decay _B+ h^l^l'^ is caused 
by the exchange of a Majorana neutrino, as illustrated 
in Fig. [U then rriip,- can be related to the Majorana 
neutrino mass m,y [5|. 

In summary, we have searched for the four lepton- 
number violating processes B'^ — >■ h~t^t^ . We find no 




(GeV/cO 



FIG. 4: 90% C.L. upper limits on the branching fraction 
as a function of the mass ?n£+^- for B"*" — >■ 7r~/j"'"/i^ (dot- 
ted/magenta line), B"*" — >■ K'jj.'^fi'^ (dash-dotted /red line), 
_B+ — )■ K'e'^e'^ (dashed/black fine), and _B+ — ^ 7r"e+e+ 
(solid/blue line). 



significant yields and place 90% C.L. upper limits on the 
branching fractions in the range (2.3— 10.7) x 10~*. The 
branching fraction upper limit for — > tt~ji~^^~^ is less 
restrictive than the result reported in Ref. |9|, while the 
— > K^fi^n'^ limit is commensurate. The limits for 
B+ A'^e+e+ and 5+ 7r~e+e+ are 30 and 70 times 
more stringent, respectively, than previous measurements 



at e~^e colliders [8|. 
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